Phantom dark energy and cosmological solutions 
without the Big Bang singularity 



A.N. Baushev a 

" Bogoliubov Laboratory of Theoretical Physics, Joint Institute for Nuclear Research 
14 1980, Joliot-Curie str. 6, Dubna, Moscow Region, Russia 



Abstract 

The hypothesis is rapidly gaining popularity that the dark energy pervading our 
universe is extra- repulsive (—p > p). The density of such a substance (usually 
called phantom energy) grows with the cosmological expansion and may become 
infinite in a finite time producing a Big Rip. In this paper we analyze the 
late stages of the universe evolution and demonstrate that the presence of the 
phantom energy in the universe is not enough in itself to produce the Big Rip. 
This singularity occurrence requires the fulfillment of some additional, rather 
strong conditions. A more probable outcome of the cosmological evolution is 
the decay of the phantom field into 'normal' matter. 

The second, more intriguing consequence of the presence of the phantom field 
is the possibility to introduce a cosmological scenario that does not contain a 
Big Bang. In the framework of this model the universe eternally expands, while 
its density and other physical parameters oscillate over a wide range, never 
reaching the Plank values. Thus, the universe evolution has no singularities at 
all. 
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1. Introduction 

Experimental data indicate convincingly an accelerating cosmological 

expansion. The substance responsible for this effect is usually called dark energy. 
If we speculate that the universe is filled with a perfect fluid with equation of 
state p — wp, the accelerating expansion appears if w < —1/3. There are a lot 
of models of the dark energy: it can be explained by a non-zero cosmological 
constant (probably, this is the most natural hypothesis, in this case w — — 1) or 
by presence of a global cosmic scalar field slowly moving down to the potential 
minimum (then —1/3 > w > —1). Until recently, the instance of w < — 1 did 
not attract considerable interest. A substance with such an equation of state (it 
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is usually called phantom energy) generally possesses physically unacceptable 
properties, for example, the sound speed in it can be super-luminal. 

Nevertheless the case of w < — 1 has currently come to the attention. First 
of all, it does not contradict to the experimental data [3, [B[ . Furthermore, it has 
turned out that small perturbations are never super-luminal at least in some 
substances with w < — 1 (see Q as an example). So, a phantom energy theory 
can be causal. 

Physical properties of a dark energy with w < — 1 differ markedly from the 
case when w > — 1. For instance, it violates the energy domination condition 
p > \p\. Since this is one of the principal assumptions of the classical black 
hole non-diminishing theorem, a black hole mass diminishes by accretion of 
the phantom energy Q. But the difference most important for the cosmology 
becomes apparent when we consider expansion of the universe containing dark 
energy. If w > —1 the dark energy density is not increasing (or even decreasing) 
as the universe expands. All bound systems (such as the solar one) remain bound 
forever; the dark energy presence manifests itself only on cosmological scales. 
If we allow w < — 1, the dark energy density grows and becomes infinite in a 
finite time. Increasing gravitational repulsion produced by the phantom energy 
will first destroy Galaxies and then any bound systems including elementary 
particles 0, This phenomenon is usually called Big Rip. However, as we will 
see below, realization of such a catastrophic scenario requires the fulfilment of 
several supplementary conditions. Even if the phantom energy prevails in the 
universe, the Big Rip does not necessarily occurs. 

There is another less evident consequence of w < — 1 : if dark energy is really 
the phantom one, the universe evolution needs not to contain a Big Bang. 



2. The phantom energy-dominated universe evolution 

The metric of a homogeneous isotropic universe can be represented as (here- 
after we use the system of units where c = H = 1, 8ttG = M~ 2 ) 

ds 2 = dt 2 - a 2 (t) [dx 2 + dy 2 + dz 2 ] = 

= a 2 (r)) [drj 2 - dx 2 - dy 2 - dz 2 } (1) 

As we have already mentioned, presence of the phantom energy in the universe 
is not enough in itself to produce the Big Rip. This phenomenon appears, for 
instance, if the coefficient averaged over the universe contents w = const < — I. 
However, it is beyond reason to believe that w remains constant during the 
cosmological evolution. In order to illustrate this let us at first consider a very 
simple situation when a universe contains only uniformly distributed phantom 
energy with an ordinary Lagrangian 




For such a field 

p = -4> 2 /2 + v (</>), v = -ftl^-v{4>) (3) 
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Hereafter a dot over a symbol denotes derivative with respect to time t. We 
obtain: 

w = = L (4) 

One can see that w is not a constant: it depends on the amplitude of the field 
and on the potential energy density V(<j>). This question has been discussed 
in close detail in Q. To summarize briefly the results: if the potential is not 
very steep (grows slower than V(<fi) oc (j) 4 ) then w tends to — f , and the density 
becomes infinite only when t — > oo. So, no Big Rip appears in this case, though 
the universe reaches the Planck density in a finite time. For steeper potentials 
a big rip singularity appears even if w — > — 1. Even the parameter w can tend 
to — oo for a v ery steep V ((/)). Finally, if V{<j>) has a maximum, w — > — 1 and 
p — > const 0, Hoj |. It is important to emphasize that a very steep potential 
is necessary to provide a constant w < — 1: for any polynomial potential, for 
instance, w tends to —1. 

As an illustration let us consider the simplest case of V(<j)) — m 2 (f> 2 /2. The 
Einstein equation can be written as 



m 2 cj) 2 



Then the equation of motion (j> + 3H(f> — V'((f>) = takes the form 

4> + 3H<j) - m 2 (t) = (6) 

Initially the phantom field rapidly grows as well as the Hubble constant H 
in accordance with ([S]). In due course, the friction term 3£f(/> becomes very 
large, and the phantom field remains almost constant. For this regime we have 
4> <C 3H<j>, 4> 2 <C m 2 2 , and (|5I6|) can be rewritten as 

^mM„^, H~—^ (7) 
V P V 3' M p ^6 V ' 

It is easy to see that no Big Rip occurs and w — > — 1. Moreover, equations ([7]) are 
closely analogous to those describing the chaotic inflation in the early universe 



In both the instances the scalar field and its energy density changes very 
slowly, and the universe evolves in a quasi-De-Sitter regime. The only difference 
is that in the case of chaotic inflation the field slowly decreases, while in the 
considering case it slowly increases (usually it is called 'slow-climb'). 

In case of another potential V((f>) the University evolution dynamics can be 
quite different, in particular, the slow-climb may be lacking. Rapid universe 
expansion, however, is a common feature of all phantom-containing models. 
Indeed, the density of the phantom matter (contrary to the 'normal' one) grows 
with time, and finally it begins to dominate. Then even in the softest case 
of w = —1 (that is not exactly a phantom energy) the scalar factor a grows 
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exponentially; if V(<p) is steep, the expansion can be much faster. So, if the 
universe contains phantom energy, it eventually comes to an inflation-like stage. 

The second effect that is able to prevent the Big Rip is a possible phantom 
energy interaction generating 'normal' particles. The normal substance gives 
positive contribution to the pressure that decreases the w value averaged over 
the universe contents. Eventually the Big Rip may be prevented. 

It is generally believed that the phantom energy does not interact with 'nor- 
mal matter', but it cannot be absolutely true, if for no other reason than the 
gravitational interaction. Let us consider a toy illustration for a start. As noted 
above, in the phantom scenario increasing gravitational repulsion produced by 
the phantom energy will finally destroy even elementary particles. We consider a 
hadron (more precisely, a quark-antiquark pair qq) ripped by the gravitational 
forces 0- When the distance between the quark and antiquark increases, the 
color field lines of force between them are pressed together into a string-like 
region. The gravitational repulsion stretches the color lines of force until the 
increasing potential energy becomes sufficient to create another qq pair. It di- 
vides the string on two strongly-coupled qq pairs, but they are also ripped up 
by the gravitational repulsion. So, this leads to intensive hadron production in 
the moments just before the Big Rip. 

Another mechanism of 'normal' matter generation in the phantom-dominated 
universe is particle production in the cosmological gravitational field. In the 
modern universe this effect is completely negligible, but if the universe is phantom- 
dominated, its density (and, consequently, the gravitational field) grows with 
time making the particle production very intensive. In [l2l ] the process was 
considered under assumption that the universe is filled with a perfect fluid with 
w = const < — 1. The influence of the 'normal' matter on the universe expan- 
sion was neglected. It was shown that if the conformal time 77 is chosen so that 
77 < 0, and the density becomes infinitive when rj — > —0, the generated 'normal' 
matter density can be written as 

Pnorm = Crf , where /3 = - - (8) 

I + 3w 

In order to examine the problem in more detail we consider a universe filled 
with a perfect fluid with w = —4/3. Initially we neglect the particle production. 
Then the cosmological equations can be easily solved. It is convenient to choose 



1 This process becomes possible when the energy of the gravitational repulsion between 
the quarks in the hadron becomes comparable with the typical strong interaction energy scale 
1 GeV. The phantom energy density at this moment depends, in particular, upon the phantom 
energy model. The gravitational repulsion force between two quarks can be estimated as 

,4ttG 

f g ~ (1 + 3w)— — pmR 

where m is the constituent quark mass, R is the typical hadron radius. Using the condition 
f g R ~ 1 GeV and taking m = 300 MeV, -R = 10 -13 cm and w ~ — 1, we obtain for the 
density of the phantom energy p ~ 10 54 g/cm 3 
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I 3 2 

the following initial conditions: pdeo, &o = \ 1 Vo — — ( we denote the 

V SnGpdeO 3 

phantom energy density by pde and its initial value by pdeo)- The value of rjo is 
taken so that the density becomes infinitive when r\ — > — 0. We obtain 

Pde carj'i (9) 

Now the particle production can be added with the help of ([8]) . We will consider 
the generated substance to be relativistic, i.e. with the equation of state p nor m = 
pnorm/3. Then the system of cosmological equations can be written as 

1 da\ 2 1 



a 2 d V J ~ Wl^ + Cr,-*) (10) 



d(p de + Cf] ?) Ida _4 

Tr, = a^- ACl1 3) (U) 

dt = adr) (12) 



The Hubble constant is given by the equation 

1 da 1 da 
a dt a 2 drj 



H ^~ = -2- (13) 



There is no accurate estimation of the C value, but we can avoid this difficulty 
choosing properly the initial density pdeo- Indeed, comparing (jHJ and (|9|), one 
can see that the density of the 'normal' matter grows faster, and eventually it 
becomes comparable with the phantom one. This is the moment of our interest. 

_ 4 

Therefore, it is reasonable to take pdeo = 100 • Ct] q 3 , i.e. to choose the initial 
moment so that the initial phantom density is only hundred times higher then 
the 'normal' one. Since the particle production is a weak effect, this value of 
Pdeo must be quite high. 

Results of the (I10j|12[) system integration with the aforecited boundary con- 
ditions are represented on Fig. |TJ|2I It is convenient to use H ■ (SiTrGpdeo/^)^ 1 ^ 2 
instead of H and t- (SnGpdeo/S) 1 ^ 2 instead oft in order to avoid the dependence 
of the graphs from the uncertain parameter pdeo- It is seen from Fig. [1] that 
the Hubble constant grows progressively slower, reaches its maximum H max , 
and begins to decrease. On the contrary, w constantly increases from the ini- 
tial value that is close to —4/3. Finally it becomes positive. Of course, both 
the effects result from the high growth fraction of the 'normal' matter in the 
universe that makes a positive contribution to the pressure. One can see from 
Fig. Q] that H m ax — 2.11 • (SnGpdeo/S) 1 ^ 2 . As we have already mentioned, pd e o 
is a very big quantity, so, H max must be also quite large. 

Extension of the solution to w = makes no sense: the equation (jHJ is 
obtained under the assumption that the dynamics of the universe expansion 
is determined only by the phantom energy. It becomes inapplicable when the 
fraction of 'normal' matter gets noticeable. It seems reasonable to contend, 
however, that in the course of the phantom-dominated universe evolution the 
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Figure 1: The time dependence of the Hubble constant H. We use H ■ (87rGp(ieo/3) 
instead of H and t ■ (SirCpdeO /3) 1//2 instead of t. 



0,0- 



-0,2- 




Figurc 2: The time dependence of w. We use t ■ (STrGpdeo/3) 1 / 2 instead of t. 

quantity of 'normal' matter becomes comparable with the phantom one, and the 
value of w grows hindering from the Big Rip. To this we can add that we have 
considered the case of w = const, which claims an extremely steep phantom 
field potential V (</)), as we could see above. For a more realistic case of the 
polynomial potential w tends to zero even in the absence of 'normal matter'. 
This makes the occurrence of the Big Rip highly unlikely. 

3. Discussion 

Of course, the suggested model is ingenuous. Intensive phantom field may 
interact in some other way decaying into 'normal' matter. Whatever the mech- 
anism of the transformation may be, the universe after it has the following 
properties: 

1. It has just passed through the stage of very rapid (at least, exponential) 
expansion. 

2. It is flat and homogeneous. 

3. It is 'normal' matter-dominated. 

4. The Hubble constant is very large. 
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It is precisely these physical conditions that existed in our universe ~ 13.7 
milliard years ago, just after the inflation. If we assume that the inflation was 
caused by the phantom field (as we could see, the equations of the phantom 
cosmology {7} are able to be closely analogous to those of the chaotic inflation 



11|), this makes possible to construct a model of universe evolution that does 



not contain a Big Bang. The universe leaves inflation being matter-dominated. 
As it expands the density of matter decreases, while the phantom density grows: 
eventually the universe passes into the phantom-dominated stage. As this takes 
place, the total density and the Hubble constant stop diminishing and begin to 
increase. Gradually the expansion becomes very fast, leading to an inflation- 
lake stage. Finally the phantom field decays into 'normal' matter, and the 
cycle repeats. Thus, the universe eternally expands, while its density and other 
physical parameters oscillate over a wide range, never reaching the Plank values. 
The universe evolution has no singularities like a Big Bang or a Big Rip. 

The foregoing model of self-reproducing cyclic universe evolution is, in a 
certain sense, a return to the Steady State theory 13, 14, 15, Indeed, in 
both the cases the universe expands forever with no Big Bang, and the matter 
is generated by the decay of some field (in our case this is the phantom field). 
The difference is that in the classical Steady State universe this process goes on 
constantly providing constant universe density, while in the considering case the 
universe evolves cyclically returning to its original state, but the density varies 
within wide limits. Presence of the evolution allows to avoid most difficulties 
inherent to the classical Steady State model. Thus, the considering model may 
be called the Steady State on the average universe. 

The presence of the inflation-like stage in the considering scenario solves 



most of the problems appearing in an inflationless cosmological theory 
The size of the universe during this stage increases so much that an observer 
in the after-inflation epoch can see only a small part of it. Consequently, the 
universe looks very flat and homogeneous. Moreover, as we could see ([7]), it is 
possible to choose the model of the phantom energy so that even the dynamics 
of the inflation-like stage is similar to the standard chaotic inflation scenario 



ill- 

Is the proposed model valid for the real universe? The answer to this question 
is currently not known with certainty, primarily because of the absence of direct 
experimental proofs of the phantom matter existence. A non-zero cosmological 
constant is quite enough to explain the experimental data. On the other hand, 
observational data show that allowed values of the parameter w of the modern 
universe (at 2a) lie between —2.25 and ~ —0.5 [5j. So, the existence of the 
phantom energy is not experimentally forbidden. 

Theoretical status of the phantom energy is also vague. There are a lot 
of Lagrangians known, which describe fields with the necessary properties. It 
is still unclear, however, whether it is possible to introduce such a Lagrangian, 
completely free from paradoxical physical effects like super-luminal sound speed, 
catastrophic quantum instability of the vacuum etc. 

Nevertheless, if the phantom energy does exist, the considering model is an 
interesting instance of cosmological solutions containing no singular points like 
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a Big Bang or a Big Rip. It is worthy of note that, though in our case the 
universe density never reaches infinity (or, more precisely, Planck density), all 
the physical parameters vary over a very wide range and amount up to extremal 
values. So, the universe passes through the stages that can be called physically 
critical. 
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